Doping of Si nanoparticles: the effect of oxidation 
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The preferred location of boron and phosphorus in oxidized free-standing Si nanoparticles was 
investigated using a first-principles density functional approach. The calculated formation energies 
indicate that P should segregate to the silicon core, whereas B is equally stable in the Si and Si02 
regions. Our models thus suggest that, in contrast with nanocrystals with H-terminated surfaces, 
the efficiency of phosphorus incorporation in oxidized Si nanoparticles can be improved by thermal 
annealing. 



Silicon nanocrystals have emerged as a promising strat- 
egy to achieve light generation and amplification us- 
ing silicon. Adding to the advantages of bulk sili- 
con, they offer variable bandgap, increased lumines- 
cence efficiency [lj], and the possibility of multiple exciton 
generation 0, Q- Together with size, shape, and surface 
functionalization, impurity doping can be used to con- 
trol the optical and electronic properties, conferring to 
the silicon nanocrystals an extraordinary degree of de- 
sign flexibility. 

Boron and phosphorus doping has been achieved 
both for self-standing nanoparticles synthesized using 
plasma approaches 0, and for Si nanoparticles em- 
bedded in glass matrix (usually Si02, borosilicate or 
phosphosilicate) @ . However, the doping efficiency is of- 
ten limited by the segregation of the impurity to the sur- 
face or to the surrounding matrix 0, [H, 0, [3 ■ 

The tendency for impurity segregation depends on the 
nanoparticle morphology and growth method. In Si 
nanocrystals with hydrogen-terminated surface, accord- 
ing to both experimental and theoretical studies, P seg- 
regates to the surface, saturating Si dangling bonds IJ, la 
[Tof . Boron was found by theoretical studies @, H, EH to 
be more stable in the sub-surface layer or at the surface 
itself, but in practice seems to be incorporated in the 
nanocrystal core, since subsequent oxidation of the outer 
shells of the nanocrystal does not decrease the fraction 
of electrically active B atoms 

The opposite behavior is observed in Si02-cmbedded 
silicon nanocrystals. In this case, phosphorus segregates 
to the Si-rich region [28| . Segregation to the silicon region 
has also been observed in Si-SiC>2 interfaces, B where it 
has been confirmed by atomistic modeling |12|. For heav- 



ily implanted material, though, P pile-up at the inter- 
face is observed instead 13[ . The case of boron is not 
so clear. It has been found to have a segregation coeffi- 
cient s = Csi/Csio 2 = 0-3, where = Cgi and Cgio 2 are 
the equilibrium concentrations @. The piling- up at the 
Si-side of the interface during segregation was justified 
by its low diffusivity in SiC>2. However, theoretical re- 
sults indicate that it is more stable at the Si-side of the 
interface or in the Si bulk.pj, Q [H[ 

In this paper, we use first principles calculations, based 
on density functional theory, to show that in Si nanocrys- 
tals covered by a SiC>2 shell, the preferred location of P is 
at the nanocrystal core, whereas B has similar energy on 
either side, thus confirming that the behavior is different 
from H-passivated silicon nanocrystals. 

The calculations were carried out using the Aimpro 
code[l§-[2lj[. The core electrons were modeled with the 
dual space separable pseudopotentials by Hartwigscn, 
Goedcckcr and Hutter [231 ■ The basis set consisted of 
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atom-centered Cartesian-Gaussian functions [23[. A con- 
tracted basis with 13 functions per atom was used for Si, 
whereas for O and B we used uncontracted basis sets with 
a total of 40 and 28 Gaussian basis functions per atom, re- 
spectively. A Pade parametrization of a functional based 
on the local density approximation was employed for the 
exchange and correlation energy [24], 

The models for the oxidized silicon nanoparticles con- 
sisted of an approximately spherical Si core of 1.5 nm 
diameter, surrounded by an amorphous Si02 shell with 
about 2 nm outer diameter. Three models were consid- 
ered. The starting structures were obtained from the re- 
sults of molecular dynamics simulations of 1800 K anneal- 
ings during 2.0 ps, 2.3 ps and 2.6 ps, as detailed in Ref . l25l. 
Dangling bonds were passivated with hydrogen atoms. 
The compositions of the Si-NPs were Sii6iOi96H5g for 
nanoparticles I and II, and SiigiOiggHsy for nanoparti- 
cle III. Nanoparticle I was the same used in a previous 
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study [H|. 

The models had no periodic boundary condi- 
tions imposed. The starting structures were relaxed at 
K, using a conjugate gradient algorithm, with the AlM- 
PRO code. The radial pair distribution function for Si, 
shown in Fig. [1] for the three nanocrystal models, shows 
that the core consisting of the first two shells of atoms 
surrounding the central Si atom has crystalline character, 
whereas the SiC>2 is amorphous. 
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FIG. 1. Radial pair distribution functions for nanoparticles I- 
III, calculated relative to the central Si atom, for the Si species 
(black lines) and O species (lighter line). Inset: Structure of 
an undoped nanoparticle (nanoparticle I). Silicon atoms in 
the core region are represented in a darker shade of gray. 

The preferred locations for substitutional boron (B s ) 
and substitutional phosphorus (P s ) were investigated by 
comparing the energy of the doped nanoparticles ob- 
tained by replacing boron or phosphorus for each of the 
161 silicon atoms in the three models I-III. We have con- 
sidered the equilibrium charge states of B s and P s , which 
are Bj and P+ both in Si and SiC^lIS]- The formation 
energies (Ef) of Bj and P+ are given relative to those 
of the equivalent defects in bulk silicon, calculated using 
as standard 512 atom supercell: 

E) = £[NP* : X] - £[NP] - {£[Si 5 ii : X] ~ £[Si 512 ]}, 

where 'NP' is the undoped silicon nanoparticle and 
NP' : X is a similar nanoparticle but with X = B~ or 
P + at the zth Si site. The energy scales of both systems 
have been aligned using the vacuum energy as standard. 

The results show an appreciable difference between B s 
and P s (Fig.[2|). For boron, the formation energy is nearly 
independent on the location. The average formation en- 
ergies arc slightly lower at the interface, but the differ- 
ence between average values is smaller than the width 
of the energy distributions represented by the errorbars. 
Phosphorus, however has a preference for the Si core re- 
gion. At the oxide shell, the formation energy is higher 
by about 2 eV. 

A fundamental factor determining the change of for- 
mation energy from the core to the shell is the differ- 
ence between the bonding energies of B or P to Si or 
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FIG. 2. Formation energy of the Si02-covered nanoparticle 
doped with B (top) or P (bottom), as a function of the dis- 
tance from the center. The filled symbols correspond to av- 
erage formation energies for distances between r — 2 A and 
r + 2 A, and their errorbars represent the statistic standard 
deviation of the samples. Open squares represent the distri- 
bution of oxygen (coincident for I, II and III). 



O. We analyze the correlation between the formation en- 
ergy and the type and number of bonds to the nearest 
neighbors, counting the number of bonds in the basis of 
a geometrical criteria: if the distance between two atoms 
exceeds the sum of their Van der Waals radii by less than 
0.5 A, we consider the two atoms bonded [27}. The top 
histogram of Fig. [3] shows that the energy distribution 
for substitutional boron is asymmetric, with a negative 
skew. In general, three-fold coordinated structures have 
low energy. Four-fold coordinated structures are evenly 
distributed on the two sides of the peak, independently 
of B~ bonding to Si or O. This shows that the cost of 
replacing a Si-Si bond by a B _ -Si bond is approximately 
the same as the cost of replacing a Si-0 bond by a B~-0 
bond. 

For P+, the energy distribution is bimodal (Fig. [3l 
bottom). The lowest energy peak is centered at about 
2 cV below the average formation energy (Ef), whilst 
the higher energy peak is about 1 eV above it. The low- 
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est energy peak is dominated by four-fold coordinated 
defects. Structures where P+ bonds to four Si neighbors 
have the lowest energy, followed by the structures where 
P + is bonded to both O and Si. The samples where P + 
has four oxygen neighbors are highest in energy amongst 
four-fold coordinated structures. Hence, the preference 
of phosphorus for the Si core is driven by the lower en- 
ergy necessary to form a P + -Si bond at the expense of a 
Si-Si bond, as compared to the energy necessary to form 
a P+-0 bond at the expense of a Si-0 bond. 
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FIG. 3. Energy vs. type and number of nearest neighbors 
(NN) for B~ (top) and P + (bottom). The formation energy 
is given relative to the average of all samples {Ef). The his- 
togram is stacked. 



We now compare the results of the present investi- 
gation to those found for H-terminated nanocrystals [8|. 
When discussing the behavior observed in samples grown 
by different methods, we assume that the dopants have 
enough thermal energy to diffuse to the most favorable 
dopant locations. During growth, that is very likely the 
case. 

Since H is less electronegative than O, in H-terminated 
Si nanocrystals the formation energy of four-fold coor- 
dinated B~ or P + is nearly independent on the lattice 
position it occupies Q. Both B and P arc stabilized at 
the surface by loosing one hydrogen to become three-fold 
coordinated, but the energy released is greater for the 
latter. This may explain why P segregates to the sur- 
face of H-terminated silicon nanocrystals, but the same 
has not been observed for bQ. On the other hand, we 
find that P should diffuse from the SiC>2 shell of SiC>2- 
covered nanoparticles or to the Si-Si02 interface to the 
Si nanoparticle core. This is in agreement with the seg- 
regation behavior observed for Si nanocrystals embedded 
in an Si02 matrix [28|. Boron, however, is predicted to 
have no energetic preference for the Si or SiC>2 regions. 
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